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^-". . We consider the addition of a large number of massive dynamical flavors to the Klebanov-Witten 

L-. theory, the quiver gauge theory describing the low energy dynamics of N c D3-branes at the conifold 

singularity. Massive flavors are introduced by means of Nf D7-branes which are holomorphically 
embedded and smeared along the transverse directions. After some general comments on the 
validity of the smearing procedure, we find the full backreacted supergravity solution corresponding 
to a particular class of massive embeddings. The solution depends on a running effective number 
of flavors, whose functional form follows from the smeared embedding. The running reflects the 
integrating in/out of massive degrees of freedom in the dual field theory as the energy scale is 
changed. We study how the dynamics of the theory depends on the flavor parameters, mainly 
focusing on the static quark-antiquark potential. As expected, we find that the dynamical flavors 
tend to screen the static color charges. 



1 Introduction 

The construction of string duals of gauge theories with dynamical flavors is a task of obvious 
interest. The gauge theories usually describe the low energy dynamics at the intersection of 
N c "color" and Nf "flavor" D-branes and evaluating the backreaction of the full system is 
not easy in general. If the flavor branes are on top of each other, the supergravity equations 
of motion cannot be simply reduced to ordinary first order differential equations in a radial 
variable. Instead, they are partial differential equations in a number of variables given by the 
number of directions transverse to both the color and the flavor branes. For D3D7 models, 
this kind of scenarios have been studied in [1] . Other localized constructions were discussed 
in 0. 

To avoid this complication one can consider simplified set-ups where the flavor branes are 
homogeneously smeared [3] in the transverse space. Whereas this restricts the possible ways 
of introducing fundamental matter (one requires that the symmetries of the unflavored theory 
should be effectively unbroken), it gives a useful framework to analyze physical features of 
the resulting flavored theories. Moreover, the "smearing trick" is not only computationally 
helpful. As we are going to argue in the following section, it is in a sense a preferred choice 
if one wants to avoid singularities and large string couplings near the position of localized 
sets of branes. 

Following this prescription, the construction of string duals of flavored supersymmetric 
theories in the Veneziano limit N c , Nf 3> 1 with Nf/N c fixed, has been possible. In jU [5J [6] 
(resp. [7]) the string dual of a SQCD-like theory with massless (resp. massive) dynamical 
flavors was found. The setup is determined by D5-branes wrapped on compact (for the 
color branes) or non compact (for the flavor branes) two-cycles. The theory is the flavored 
version of the confining Chamseddine-Volkov-Maldacena-Nunez (CVMN) solution [8]. In [9] 
a large number of massless flavors was added to the conifold conformal theory of Klebanov 
and Witten (KW) [10J. In [11] the construction was extended to the confining Klebanov- 
Strassler [T2] case. Other related setups and studies can be found in [T3 l 1^ 1 fT5l fT6] . 

In this paper we focus on the flavored Klebanov- Witten model and extend the analysis of 
[9] by considering the case where the dynamical flavors are massive and all with the same 
constituent mass m q . Fundamental flavor multiplets are added to the theory by means of 
D7-branes, which are wrapped on non-compact 4-cycles and holomorphically embedded in 
the background in order to preserve the M = 1 supersymmetry [17]. The authors of [9] 
found the supergravity solution generated by a "massless", smeared D3D7 system. This has 
a running dilaton and non trivial F\ (sourced by the D7-branes) and F 5 RR fluxes. The near 
horizon geometry has no AdS factor and this translates in a broken conformal invariance in 
the dual field theory. The running dilaton blows up at a certain value of the radial coordinate 
corresponding to a field theory UV Landau pole. Moreover, the background has a curvature 
singularity when the dual gauge theory is in the far IR. 



In the following sections we will show how with massive flavors one can find a IR regular 
solution. Of course, the latter will not alter the UV behavior of the field theory and the 
Landau pole will still be present. The flavor D7-branes are embedded in such a way that they 
reach a minimal radial distance p q from the bottom of the space. This distance is related 
to the mass parameter m g , just as the field theory energy scale is related to the radius p. 
Energies larger (smaller) than m q map to radii larger (smaller) than p q . The knowledge of 
the conifold geometry allows to find the density distribution of the smeared flavor branes as 
a function of the radial coordinate. 

With the density distribution explicitly calculated, we are able to solve the supergravity 
equations of motion coupled to the Nf D7-brane sources. The net result is interesting. As 
it was anticipated in [9], the first order equations for the background fields following from 
the supersymmetric fermionic variations retain the same form as those in the massless case, 
modulo a substitution of Nf with a running effective number of flavors, Nf(p). The function 
Nf(p) is related to the flavor density distribution mentioned above. This effective running of 
the number of flavors has a nice field theory interpretation. As we go towards energy scales 
much larger than the mass m q the theory resembles the massless one; at energies lower than 
m q , the flavors can be integrated out and the theory looks like the unflavored one. We will 
see that the function Nf(p), as well as the whole supergravity solution we will deduce from 
it, precisely reproduces this field theory expectation. The function Nf(p) in fact has the 
shape of a smoothed out Heaviside step function NfQ(p — p q ). In contexts where extracting 
the function Nf(p) from the smeared massive embedding could be technically difficult, it 
is useful to plug in the supergravity equations just the simple Heaviside function. In this 
paper we will compare the results obtained by using the correct Nf(p) with those deduced 
using the step approximation, finding that, in the flavored KW model at hand, the latter 
works quite well, at least at small m q . This approximation was used to study novel effects 
of massive dynamical flavors in the confining C VMN background in [7] . 

Once the full flavored KW background has been found, one can use it to study how the 
strongly coupled dynamics of the dual gauge theory is affected by the dynamical flavors. One 
of the expected effects is the screening of the color charge. In this paper we use standard 
holographic techniques to extract the static potential between an external (i.e. extremely 
massive) quark- ant iquark pair Q, Q probing our flavored gauge theory. The potential results 
to have a Coulomb-like shape and its behavior as a function of the sea quark parameters 
Nf,m q precisely accounts for the expected screening of the color charges. We also study 
the behavior of the screening length as well as that of the minimal quark-antiquark distance 
at which the QQ bound state can decay into a pair of specific heavy-light bound states 
Qq + qQ (by popping out from the vacuum a dynamical quark-antiquark pair q, q). Our 
analysis shows that the screening length is an increasing (resp. decreasing) function of m q 
(resp. Nf). The above mentioned minimal distance, which we call "string breaking distance" 
L s f), is instead a decreasing function of both Nf and m q . 



The paper is organized as follows. We start in Section [2] by reviewing the smearing tech- 
nique, pointing out which are its limits of validity and its relevance. In Section [3] we consider 
certain generalized massless and massive D7 embeddings on the conifold. Their smearing is 
considered in Section H] where the expression for their radial distribution density in a massive 
case is calculated. In Section [5] we find the supergravity background dual to the KW theory 
coupled to massive dynamical quarks and discuss the regimes where we can trust the solu- 
tion. In Section O we rewrite the supergravity solution by using the Heaviside step function 
approximation. We then study, in Section [TJ the dependence on the flavor parameters of the 
static quark-antiquark potential and of the screening and string breaking lengths, making a 
comparison between the results found using the "true" supergravity solution and those fol- 
lowing from the step function approximation. In Appendix |A] we review the orbifold origin 
of some D7 embeddings on the conifold. In Appendix [B] we study the dependence of the 
static quark-antiquark potential on a particular integration constant. 

2 Comments on the smearing technique 

The technical trick we adopt to construct the string duals of gauge theories with a large 
number of dynamical flavors, consists in considering the flavor branes as homogeneously 
smeared along the transverse directions [3j. This trick allows one to find relatively simple 
solutions having the largest possible degree of symmetry and taking into account the full 
backreaction due to iV c "color" and Nf "flavor" D-branes. Moreover, as we are going to 
argue, the trick is in a sense forced by the approximations used to find the full solution. In 
fact, the starting point on the low energy string side is an action of the form [18] 



N f 

S = S II + J2[S l Bi + S iWZ ] J (2.1) 

i=l 

where the first term is the (IIB or IIA) bulk supergravity action and the remaining ones 
are the Dirac-Born-Infeld and Wess-Zumino part of the flavor brane action, with the flavor 
branes embedded in the background. One usually first solves for the embedding equations 
and then for the remaining bulk equations on-shell. 

Let us now consider Nf coincident flavor branes, such that the corresponding flavor sym- 
metry group is U(Nf). In this case one expects corrections to the DBI action coming from 
the fact that an open string can end on the Nf branes, so that the effective coupling is g s Nf. 
In the quenched approximation Nf <C N c and the coupling is small, but if Nf = 0(N C ) 
the string theory would be strongly coupled. The DBI action was shown to appropriately 
describe the g s iV/-leading order dynamics of open strings in a generic background [IH] . Thus, 
if g s Nf is order one, the corrections to the DBI can be large. This fact has a precise analogue 
in field theory in the Veneziano limit [20]. At any given order in N c , the insertion of n quark 
loops ("windows") is governed by the parameter {gyM^f) n an d can be done perturbatively 



only for small g^M^f- If the latter is not small, the sum over any number of windows must be 
performed non-perturbatively. In the brane language, the "one window graph" corresponds 
to the DBI contribution, which is then a good approximation only for small g s Nf. 

Crucially, in the smeared setup the effective coupling g s Nf is further suppressed^ Due to 
the smearing, the flavor symmetry U(Nf) is generically broken to U(l) Nf and the distance 
between two generic flavor branes is large (in string units) in the supergravity approximation. 
Now, a typical string process will involve a space-time region of size of order one in \fa' units: 
in such a region, thus, only a small fraction of the Nf smeared flavor branes will be available 
for the process. If we denote as R the typical radius of an internal dimension of the geometry 
in y/cV units, the number of flavor branes involved in a typical process will be of order Nf/R d , 
where d is the codimension of the flavor branes (the number of dimensions involved in the 
smearing). Since in the supergravity approximation R ^> 1, the effective coupling to the 
flavor branes (g s Nf)/R d will typically be small even if Nf = 0(N C ). Thus, in cases as jl]- 
[7j , [13] , [15] where the supergravity regime can be attained with no restrictions on Nf/N c , 
the use of the DBI action is still justified also for g s Nf of order one. For the case considered 
in this paper, nevertheless, we will show in section 15.51 that the validity of the supergravity 
approximation implies Nf <C N c and that the validity of the DBI does not impose any further 
restriction. 

A more detailed discussion of these issues will be presented in [2T1 . 



3 D7 embeddings on the conifold 

The low energy dynamics of N c D3-branes at the conifold singularity 

z x z 2 = z 3 z 4 , (3.1) 

where the z-i are complex coordinates, is described by an M = 1 superconformal quiver 
gauge theory with gauge group SU(N C ) x SU(N C ) and bifundamental matter fields Ai,A 2 
and Bi,B 2 transforming respectively in the (N C ,N C ) and in the (N C ,N C ) representations of 
the gauge group [10]. The matter fields form two SU(2) doublets and interact through a 
quartic superpotential 

W KW = e ij e kl [A i B k A j B l }. (3.2) 

Here and in the following, traces over color indices are implied. In the N c = 1 case it is 
not difficult to show that the moduli space of the theory is in fact described by a conifold 
singularity, with the following map between geometrical data and mesonic vevs 

z 1 = A 1 B 1 , z 2 = A 2 B 2 , z 3 = A x B 2 , z 4 = A 2 B X . (3.3) 



1 We are grateful to Carlos Nunez for this crucial observation and for his relevant contributions to the 
discussion presented in this Section. 



The conifold is a 6d Calabi-Yau cone over the T 1,1 Sasaki-Einstein manifold. Its Ricci flat 
metric is usually written as 

ds 2 c = dr 2 + r 2 ds^i,i , 

d4u = \JZ[ de2 i+^ 2e M] + \W + iZcose t d^]\ (3.4) 

where the range of the angles is ^ G [0, 47r), <^i £ [0, 27r), 6>j G [0, it]. In terms of these 
coordinates we can write 

= r f e | W-vi-va) sin ^1 s in - , Z2 = r I e l GHiPi-H*) cos ?i cos ^ 

2 2' 22' 

2 3 = r f e *(iHiw-¥») cos ^1 s in - , 24 = r le^-^ +V2) sin - cos — . (3.5) 

For our purposes of generalizing the background to an unquenched setup, it is more conve- 
nient to work with a redefined radial coordinate 

r = r e p , (3.6) 

where p is a dimensionless quantity and r is a constant which, in order to simplify notation, 
we set to one from now on. The supergravity background sourced by the N c D3 branes has 
a metric of the known warped form 

ds 2 = /i- 1/2 (p)[cfe^] + h 1/2 (p)[e 2p dp 2 + e 2 'd4i,i] , (3.7) 

constant dilaton <ft, and N c units of F§ RR flux through T 1 ' 1 . In the decoupling limit, where 
the geometry is dual to the Klebanov-Witten fixed point, the relevant metric is AdS§ x T 1,1 . 

The D7-branes used to flavor the KW theory, are taken to be extended along the Minkowski 
4d directions and along a non compact four- dimensional submanifold in the transverse space. 
This makes relatively easy the task of solving for the embedding equations. In fact the warp 
factor drops out in the D7 DBI action and the spacetime effectively seen by the brane is a 
direct product of Minkowski spacetime and the conifold. This allows to write the embedding 
equations in a relatively simple way in terms of the conifold coordinates Zj. A detailed study 
of the holomorphic D7 embeddings on the conifold can be found in [221 123| EH [25] . 

Two classes of embeddings, having z\ = \i and z\ — z^ = \x as representative elements, 
were mainly considered in the literature. On the field theory side the two classes correspond 
to the addition, to the KW superpotential, of extra cubic and quartic terms, respectively, as 
well as of mass terms for the fundamental flavors when /i^O. In the following we will focus 
on the first class only. In Appendix [A] we will briefly review the field theory interpretation 
of both kind of embeddings, starting from an orbifold construction. 



3.1 Massless embeddings 

The simplest embedding we want to consider is z\ = 0. We call this embedding massless 
since it is extended down to p = — oo. From ( 13. 5ft above, it obviously implies 9\ = or 
9 2 = 0. Thus, there are two branches. In each branch, the D7 brane fully wraps one of the 
two two-spheres of T 1,1 and chooses a point in the other one. Notice that usually choosing a 
point in a S 2 requires two equations (fixing 9, <p), but of course if one is at a pole (sin 9 = 0), 
it is not necessary to specify (p. Thus we can describe a massless embedding of this kind by 
choosing two points, one for each sphere, fixing the position of each branch. 

From the field theory point of view the two branches correspond to two classes of flavor 
multiplets, which we denote by q x , q\ and q 2 , q 2 . The addition to the Klebanov-Witten theory 
of Nj D7-branes with z\ = embedding, modifies the field theory superpotential as [22] 

W Z1=0 = W KW + h x q x A x q 2 + h 2 q 2 B x q x , (3.8) 

where, here and in the following, sums over the Nf flavor indices are implied. The classical 
flavor symmetry preserved by W Zl =o is U(Nf) x U(Nf). Because of the SU{2) x SU{2) 
symmetry, one should be able to pick a (two-branch) embedding by picking any point on 
each of the spheres. This means that we should have a four-parameter family of massless 
embeddings. This family is described by the more general embedding equation [22] 

4 
E"A = 0- (3-9) 

i=l 

with the complex constants a, spanning a conifold (up to overall complex rescalings) 

a x a 2 — a^a 4 = . (3.10) 

Notice that embeddings like Z\ — z 2 = are not in this familyo In (13.91) . (13.10J) . we can 
rescale the coefficients to fbo a x = 1. Then (13.101) fixes a 2 - Thus, this class of embeddings 
is parameterized by two complex numbers a%, a^, as we expect in order to be able to choose 
a point in each of the spheres for each branch. Indeed, the embedding equation 

z x + a 3 a 4 z 2 + a 3 z 3 + a 4 z 4 = , (3.11) 

nicely factorizes into two branches 

(sin y + a 3 e iipi cos j J (sin j + a 4 e iip2 cos y J = , (3.12) 



2 A SU(2) x SU(2) rotation, in fact, maps z\ — Z2 = into a generalized embedding equation as in (3.9), 
with the parameters spanning a unit 3-sphere: o.\ = —0.2, a^ = 0.4, \ai\ 2 + | CK3 1 2 = 1. 

3 Obviously this does not work for the a\ = case, which implies that one of the branches is at the south 
pole of one sphere. This is a zero measure subset of the embeddings considered. 
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so, as expected, a^, on determine at which point each branch is sitting in each sphere: 

0\ = 0^ = 2 arctan |a 3 | , ipi = <pf = n — arg[a 3 ] , 

9 2 = 9™ = 2 arctan |a 4 | , tp 2 = ^2° = n ~ ar g[ Q; 4] • (3.13) 

The constants 8f°, iff denote the position of each branch as p — ► oo (a notation that will be 
useful in the next section). We can rewrite (13. 12)) as 

Tir 2 = 0, (3.14) 

where we have defined 

ooo n noo a 

Tt = cos ^- sin - - j^-vt) sin Zi_ cos Zl (* = 1,2). (3.15) 

2 2 2 2 v ; v ; 

For later reference, it is worth remarking that 

0<|r*|<l. (3.16) 

3.2 Massive embeddings 

We now want to give a mass to the flavors. This deforms the embeddings in such a way that 
the two branches of the massless embeddings merge (the U(Nf) x U(Nf) flavor symmetry 
is explicitly broken down to U(Nf)). Let us again start with the simplest embedding z\ = 
6 2pi e l @, where p q , /3 are real numbers. This is referred to as the unit winding (m = n 2 = 1) 
embedding in the notations of [24j and it is explicitly given by 

^ = ^1 + ^2 + 2/3, e 3 ^= e / 9 (3.17) 

(sin |) (sin |) 

Clearly, p q is the minimal value of p reached by the brane. This is a connected embedding 
which for large p goes to the two branches of the corresponding massless embedding. 

We now want to consider the same generalized class for the massive embeddings as we 
used in section 13.11 The defining equation is [22] 

Z\ + olzol^Z2 + 0:3-23 + OJ4Z4 = const . (3.18) 

In terms of the angles, we can write it as 

e §P e *(iMw-va) ri r 2 = e^ p " e i(3 , (3.19) 

with the same T's defined in (13. 15)) . Something obvious that is worth noting is that the 
embedding equation is not independent on ip any more, as expected because a mass term 
breaks the U(1)r symmetry. Shifting ip by a constant shifts the phase (3. Consistently, there 



is a I in front of ip such that a 4n shift which takes us to the same point, shifts the mass 
phase by 27r, thus leaving it invariant. 

Thus, the embedding of a brane within this family depends on 6 real parameters: 9f°, ip™, 
0™, <^£°, ft, p q . In the field theory, the first four should be identified with the couplings in 
the superpotential while the last two should be related to the phase and the modulus of the 
mass term. 

The generalized embeddings we have found above, can be mapped to a field theory super- 
potential of the form 

W = W K w + hiqi[Ai + aiA 2 ]q 2 + h 2 q 2 [B\ + a 3 B 2 }qx + m 1 q 1 q 1 + m 2 q 2 q 2 , (3.20) 

where we put generic mx,rri2 mass terms (then we will want mx = m 2 ). We can formally 
rewrite W in a compact form as W = Wkw + q~Mq where M is the effective mass matrix 
for the flavors. Following the arguments in [22] , the equation det M = should correspond 
to the D7 embedding equation. Indeed, it reads 

hxh 2 [Ax + a 4 A 2 ][Bx + a 3 B 2 ] = m x m 2 , (3.21) 

which is equivalent (using, for the iV c = 1 case, the maps (J3.3P ) to the equation 

m 1 m 2 
Zx + a 3 a 4 z 2 + a 3 z 3 + a 4 z 4 = -— — . (3.22) 

hxh 2 

This is exactly of the form we proposed above, for the massless mx = m 2 = (13. lip and the 
massive (13.181) case. 

4 Smearing the embeddings 

We want now to consider Nf ^> 1 branes suitably distributed within the family of embeddings 
(I3.19p . In particular, we want to restore effectively the SU(2) x SU(2) x U(l) symmetry of 
T 1 ' 1 which is broken by flavor branes on top of each other. 

Accordingly, we want to place branes at different values of O™, ff, Of, <pf,ft, while, in 
principle, we are interested in keeping p q fixed (clearly, it is also possible to implement a 
distribution of Pq's). The symmetry preserving density of branes is 

Ofloo^ooaoo^ooa = = — sin 9?° sin 9%° . (4.11 

Pe x Vl e 2 <p 2 p d9Tdyrd9?dy¥dft 2vr (4tt) 2 x 2 ' l ' 

which is normalized as 

r de? r d^ r d9? r ^~ r ^ P9ri/Tervrf} = N f . (4.2) 



The purpose of this section is to compute how the charge density produced by this D7-brane 
distribution modifies the Bianchi identity for F\. This coupling comes from the WZ term, 
which can be written as (cfr. eq (3.1) of [9]) 



N JMs JMiq 



S WZ = T 7 '£I C 8 -+T 7 I OAC 8 , (4.3) 

yielding: 



dF 1 = -T 7 (2n 2 {w) )n = -g s Q . (4.4) 

4.1 General procedure 

Let us start by considering a generic distribution of D7-branes. The Q defined above is built 
from the orthogonal planes to the submanifolds where the branes are sitting. Each brane is 
described by two equations 

f 1 (p l]X n = 0, / 2 (p i ;x A4 ) = 0, (4.5) 

where the pi are some parameters and x M are the spacetime coordinates (// = 0, ..., 9). Locally, 
the orthogonal plane to a single embedding is described by a two-form: 

5(fi)S(f 2 )df 1 Adf 2 . (4.6) 

We now need to sum over all the branes. This becomes just an integral over the parameters 
which define the embeddings, suitably weighted with the density: 

n = Jp Pi (p l )(5(fMf2)df 1 Adf 2 )d Pl . (4.7) 

4.1.1 A simple example 

To see how the calculation is performed, let us work out the expression (14.71) for a case in 
which we already know the answer: the smeared massless embedding in the KW case, as 
discussed in [9]. There, since each embedding is disconnected into two branches, we have to 
compute separately each contribution. 

Consider the collection of branches described by 

h(9?, cpT; a?) =9i-ff?, f 2 (6?, <p?; xT) = <pt - <p? . (4.8) 

The distribution density is 

p 0rvr = -^ML9? . (4.9) 

We plug this in ([4.7]) to get 

Q = ( [ ^L sin e^°5(9 1 - 9™)5{ Vl - cp™)d9™d(p™) d9 1 A dcpt = ^ sin 9 1 d9 1 A d Vl . (4.10) 

\J 4:71 ) 4:71 

From the other branch, we get a similar contribution with (1 — > 2), so Q = -^(sin9id9i A 
dif\ + sin9 2 d9 2 A d(p 2 ). This, together with (14.41) above, agrees with equation (2.13) of [9]. 
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4.2 The Bianchi identity for the smeared massive case 

We now want to compute the expression (14.71) where /i, f 2 can be read from (13.191) and the 
density is written in ( 14.11) . Let us split the modulus and the phase in ( 13.191) . From equating 
the phases, we get 

h=il)-<fx-<P2 + 2 arg(r x ) + 2 arg(r 2 ) - 2/5 - Ann , (4.11) 

where n is any integer number. 

From the modulus (squared) we find 

/ 2 = e 3p |ri| 2 |r 2 | 2 -e 3p9 , (4.12) 

where 

IT^I 2 = - [1 - cos0°° cos^, - sinfli sinflf cos(^ - pf)] . (4.13) 

We now have to insert these expressions, together with ( 14.11) . into ( 14.71) . Let us start by 
making the integral in /3. Since the only dependence from /3 in the integrand comes from 
S(fi), we have / 5(fi)d{3 = | (notice that /i changes by Att when (5 goes from to 2ti, so, for 
any values oiip,9i,(fi, 9°°, ip°°, the argument of the delta function is zero exactly one time in 
the j3 range of integration). Thus: 

Q = —^ f (sin0~ sin9 2 x >6{f 2 )df 1 A df 2 ) d9™d V ™d9™dp 2 X) . (4.14) 

(47r) d J 

Notice that, because of (I3.16p . f 2 can only vanish if p > p q . Therefore, Q(p < p q ) = 0, as 
expected, since the flavor branes are extended in p > p q only. 

The integral (14.141) seems extremely difficult to compute. However, the result must be 
fairly simple, just of the form displayed in (3.99) of [9] for a certairo Nf(p) (a dot will denote 
derivative with respect to p throughout the paper): 

Nf(p) Nr(p) 

ft = J v ' (sin 6id0i A dip! + sin 9 2 d9 2 Adif2) J -^dp A( dip + cos 9 1 dip 1 + cos 9 2 dip 2 ) . (4.15) 

An Atc 

This is the only possibility for an exact two-form preserving SU(2) x SU(2) x U(l)^ x Z 2 

(the Z2 interchanging the 1 and 2 S 2 's). This observation will allow us to obtain an exact 

simple expression for Q since this task is just reduced to computing a single function Nf(p). 

Let us expand Q = ^flMNdx M f\dx N , where the M, N are the coordinates p, ip, 9%, 9 2 , <px, (p 2 . 

We can easily isolate the different terms of this expression and compute them in turn. The 

simplest component is 

n^ = -j^ J (sm9?sm9?8(f 2 )3e 3p \T l \ 2 \T 2 \ 2 ) d6? d<p? d9? dy? . (4.16) 



4 Notice that Nf will denote the fixed number of flavor branes while Nf(p) is a non-trivial function, 
representing the effective number of massless flavors at a given energy scale. These two quantities should 
not be confused, they are only equal in the massless case, p q — —00. 
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Figure 1: The function Nf(p) for three different values of the mass parameter (resp. from 
left to right) p q = — 5 , — 4 , —3. 

Still, this is not so easy to compute due to the delta function. But, since fl p ^ cannot 
depend on 0i,<fi (see (14.151) ). we can fix 6\ = 9 2 = such that, from (14.131) . we have 
|rj| 2 = |(1 — cos Q°°). Changing variables to a, = |(1 — cos 9°°) (notice that < a^ < 1), the 
integral above now takes a simple form: 



n 



pi/> 



TV r 
-~r~ 5(e 3p a 1 a 2 - e 3pq )3e 3p a ± a 2 da x da 2 
An J 



4ir 
Comparing the result to (14.151) we get the sought result: 



3iV> 


P 3p 9 - 


-3p f 1 
Je 3 t>q- 




da 2 


47T 




-3p 


a 2 



(4.17) 



N f (p) 

N f (p) 



3N f e 3pq - 3p (3p - 3p, 



QJ ) 



iVf 



1 - e 3p "- 3p (l + 3p - 3p,j 



(P > Pq) 



(4.18) 



whereas Nf(p) = at p < p q . Consistently, iV/(oo) = iV/ and Nf(p q ) = 0, so Nf(p) is a 
continuous function which asymptotes to the one of the massless case when p ^> p q . Notice 
also that Nf(p),Nf(p) > as required on intuitive grounds and also by the supergravity 
equations of motion pj. As it is shown in figured! Nf(p) has the shape of a smoothed-out 
Heaviside step function. 

This completes the computation of the Q since the symmetry constrains the rest of the 
components to be given by (14.151) . We have checked the result (14.181) by calculating explicitly, 
by numerical integration, the other components. 

The effective running of the number of flavors has a natural field theory interpretation. 
At energy scales larger than the mass of the matter multiplets, the theory looks like the 
massless one. At energies lower than the mass scale, the massive flavors can be integrated 
out and the theory resembles the unflavored (Nf = 0) one. At scales around p ~ p q , the 
non-trivial profile of the function should take into account threshold effects. 
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5 The backreacted solution with massive flavors 

Knowing the precise expression for Nf(p), we can find the backreacted solution dual to the 
massive flavored KW theory using the expressions of [9J. The ansatz for the Einstein frame 
metric 



ds 2 



i i / e 2g e 2f 

h *dx 2 3 + h* e 2f dp 2 + — ^ ( d °i + sin 2 9id(p 2 ) + — (dtp + J2 cosM^i) 

\ " i=l,2 " i=l,2 



(5.1) 

is a deformation, driven by the functions f(p),g(p), of the warped conifold considered up to 
now (for which f(p) = g(p) = p). The radial coordinate p varies in (— oo, pi] with the lower 
(resp. upper) extremum being mapped, via the holographic radius/energy relation, to the 
extreme IR (resp. UV Landau pole) of the dual field theory. The ansatz for the RR-forms 
and the dilaton reads 

<t> = <P(p) , 
F 5 = 277rg s N c a' 2 e-' i ^ f h{p)- 5 ^[e x0xlx2x3p -e ewe2 ^), 

Fi = s (dip + cos 9\ d(fi + cos 8 2 d(p 2 ) , 

4:71 

dF 1 = - 9s /^ (sin 6i d0! A dip! + sin 6 2 d$ 2 A d<p 2 ) + 

4:71 

1 9s s ^ P 'dp A (dip + cosOidipi + cose 2 dip 2 ) ■ (5.2) 



4tt 
The vielbein is given by 

e x ' = h- 1/4 dx\ e p = h 1/4 e f dp , 

e e ' = -^h 1/4 e g d6 t , e** = ^=h 1/A e? sm9 l d^ i , 

e i> = h 1/4 e f (dip + cos 6!d(pi + cos 6 2 dip 2 ) . (5.3) 

It is crucial to notice that the projections that define the Killing spinors have the same 
form for the deformed or undeformed ansatze [9]. Then, it is not difficult to check that 
the ^-symmetry analysis of [21] is easily generalized to the backreacted solution, and (13.171) 
is obtained without modification (when one uses the p-coordinate). This means that the 
computation of section H] and in particular the result (14.181) is directly applicable to the 
backreacted case. 

The other crucial observation is that passing from Nf to Nf(p) does not alter the form 
of the first order equations for the functions appearing in the ansatz. This is due to the 
fact that the supersymmetric fermionic variations only contain the forms F and not the dF 
terms. Thus, although the Bianchi identity for F\ is modified in this massive setup w.r.t. 
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the massless one, we can just use the first order equations found in [9] with Nf — > Nf(p). 
We thus have: 



9 = e 2/ - 2s , 
/ = 3-2e 2f ' 29 
■ _ 3g.Nf(p) * 



3g s N f (p)^ 



e T 



•>7T 



4vr 
h = -27ng s N c a' 2 e~ 49 . (5.4) 

These equations have to be solved separately for p < p q and p > p q . Then, one has to 
demand continuity at p = p q . This will be the content of the upcoming subsections. 

The constituent mass m q for the dynamical quarks can be related to p q by defining it, as 
usual, as the energy of a straight string stretched along the radial direction from p q to the 
bottom of the space, that is at p = — oo in our case: 






5.1 Region 1: p < p q 

In this region, one just has the well known unflavored system, after inserting Nf(p) = in 
(15. 4p . The general solution for 0, g, f can be given in a simple form: 

e 4> = e 4>, R 



oil 



& = c 3 (e 6p + ci 
e / = c 3 e 3p (e 6p + c x ] 
h(p) = 27n g s N c a' 2 (c 2 + ["" e" 49 ^^^) . (5.6) 



For the integration constants, we have used a notation similar to [9]. 

5.2 Region 2: p > p q 

The relevant system of equations for this region comes from inserting Nf(p) as given in (14.181) 
in (15.41) . Remarkably, it turns out that the functions f,g,<p can be explicitly integrated in 
this region too. Just as in the massless case, the dilaton is running and blows up at a 
certain p^. Fixing one of the constants of integration, which just amounts to shifting p, p q , 
we require, as in [9], that the dilaton divergence (i.e. the field theory UV Landau pole) is 
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located at pi = 0. The solution is: 

47T 



e / = c 3 



g,N f (3p - e 3 ^-3p(3/> 9 - 3p - 2) + e 3 ^(3p g - 2)) ' 

(-6e 6p p + 2e 3p+3p *(3p q - 3p - 2) - 2e 6 ^+ 3 ^ (3p g - 2))2 



1 5 
|3 



{e 6p (l - 6p) + 4e 3 P+ 3 ^(3p g - 3p - 1) - 2e 6 ^+ 3 ^(3p g - 2) + c x ) 
e 9 = c 3 (e 6p (l - dp) + 4e 3p+3p " (3p q - 3p - 1) - 2e 6p+3p * (3p g - 2) + c^ 

h{p) = 27rtg s N c a' 2 (c 2 + Te-^^dp*) . (5.7) 

5.3 Patching the solutions at p = p q 

Out of the seven integration constants appearing in (15. 6p . (15.71) . four are fixed by demanding 
continuity of the functions at p = p q , namely: 

47T 



^ s iV / ( e 3 P 9 (2-3p g )-(2 + 3p g )) ' 
Cl = - e 6 ^ +g 1 (2e 3 ^(2-3p g )-2(2 + 3p q 



c 3 = c 3 (2e 3 ^(2-3p, 2 )-2(2 + 3p (2 



~C2 = c 2 + [°e- 49 ^dp* . (5.8) 

Notice that the infrared dilaton blows up, signaling the breaking of the validity of the su- 
pergravity approximation, for p q — > 0, i.e. for dynamical masses close to the Landau pole. 

5.4 A prescription for the integration constants 

Apart from depending on the parameters N c , Nf, p q , which have a clear field theory interpre- 
tation (p g being related to the modulus of the quark masses), the solution described above 
also contains three independent integration constants 5i, c 2 , c 3 . In this subsection we will 
discuss their physical meaning and give a prescription to fix them when comparing different 
solutions, as we will do in section [71 

Let us start with c\, which severely affects the IR behaviors. If c.\ = one recovers in 
the IR limit the standard regular solution for D3-branes on the conifold. If Ci > the 
solution is of the form discussed in [26] : in this case the transverse space to the D3-branes 
is a deformation of the (Z 2 orbifold of the) conifold where a 4-cycle has blown up. The D3's 
are actually smeared along that cycle and their backreaction gives rise to a singular lOd 
background. Here, we will not be interested in such solutions and thus, by demanding IR 
regularity, we impose 

ci = , (5.9) 
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which, from (j5.8p . sets 

d = -e 6p9 . (5.10) 

The constant c 2 affects mainly the UV and corresponds to turning on a source term for an 
irrelevant operator. We will require our solutions to share a common UV behavior by setting 

c 2 = 0, (5.11) 

such that h precisely vanishes at the Landau pole. We then have <p\ p= o = 00, h\ p=0 = 
e f \ p=0 = 0. 

Finally, C3, which simultaneously rescales e 29 , e 2 f , h~i, can be reabsorbed by a rescaling 
of the Minkowski coordinates. We will adopt an ad hoc prescription and, as before, require 
that all solutions behave similarly in the UV by demanding e 9 \ p= o = 1. This fixes: 

C3 = (1 - e Qp " + 6p q e 3p ")-^ . (5.12) 

Inserting this value in (15.81) we obtain: 



'2e 3 Pi(2-3p q )-2(2 + 3p qJ 



~ C3 ~[ l-e^+6/v^ J • (5 - 13) 

Let us close this section by inserting these prescriptions in the definition of the constituent 
quark mass (j5.5p : 



in 



-jP'I 



m a = J , r . (5.14) 

2na'\ g s N f (1 - e 6 ^ + 6p g e 3 ^)e(2 e 3 ^(2 - 3p q ) - 2(2 + 3p g ))3 

Notice that lim Pq ^om q = 00. Of course we cannot trust the supergravity formulas in this 
limit since the dilaton is blowing up. The best we can say is that these results suggest that, 
within the choice of integration constants adopted above, we can explore a wide range of 
dynamical mass parameters, from m q — > to some higher m q < A^y. 

5.5 The Ricci curvature and regime of validity 

The validity of the supergravity approximation that we are using, requires the curvature 
invariants of the string frame version of (15. ip to be small in a' units. Let us focus, in 
particular, on the Ricci scalar. Using the BPS equations (15.41) . one obtains 

R S = - 2 39 '"'W h -h-*+t (l + ^'"M e»-*f+* + l e ^-V^M\ . (5.15) 

47r \ 4tt 4 N f (p)J 

Obviously, it vanishes in region 1 (p < p q ) where the near-horizon (p — > —00) metric is 
AdSs x T 1 ' 1 and the dilaton is constant. Validity of supergravity in region 1 requires Xir = 
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Figure 2: The curvature scalar as a function of p for, from left to right, 
p q = -20,-15,-10,-5,-3. 

g s N c e^ IR to be large and g s e^ IR to be small, as usual. From (J5.8P we see that these conditions 
are satisfied, at fixed p q , if N c ^> Nf ~^> 1. 

Let us now concentrate on region 2. Since e^ ~ (g s iV/) _1 and /i ~ g s N c a' 2 , it is easy to see 
that one can write Rs as J Nf/N c a'~ l times a factor which only depends on p, p q (once the 
integration constants c\, c<i have been fixed, notice that Rs does not depend on C3). Figure 
|2] depicts some examples, where the prescription of section [53] for ci, C2 has been used. 

In region 2, in order for the supergravity approximation to be valid in a large range of 
p, we need Nf <C iV c so that the Ricci scalar stays parametrically small except, eventually, 
very near the Landau polejj Notice also that the only singularity of the solution is located 
at the Landau pole (a pathological behavior expected on physical grounds). As stressed in 
the introduction, the smearing of massive flavors allows to smooth out IR singularities. 

Apart from checking the validity of the supergravity approximation which requires Nf <C 
N c and g s -C 1, we should also check the validity of the smeared approximation and the DBI 
action, along the lines of Section [2J The branes we are smearing have codimension d = 2 and, 
as we have seen, the typical length of the transverse space is governed by R ~ (N c /Nf)^. 

Within an area of a' size there are Nf/R? « Nf N c 2 flavor branes, a number which should 
be larger than 1 for the smearing to be a good approximation down to a' scales. On the 

other hand, for the DBI to be valid, g s N? N c 2 <^C 1. If we take g s ~ N~ x or Nj , it is clear 
that this condition does not impose any further restriction. We thus get the following range 
of parameters: 

K N? < N f < N c . (5.16) 

Similar restrictions on the parameters should be considered also in the the massless case 
p q = — 00. It should be interesting to investigate if there are regions of p for which our 
solution can be valid for more general values of Nf/N c and to check how the curvature - 



5 This restriction is typical in D3D7 systems pQ. Notice that, despite the limit Nf <§; N c , we are not in 
the quenched approximation, since the backreaction of the flavor branes is taken into account. 
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and, so, the validity range of the supergravity approximation - is affected by the choice of 
the integration constants. 

5.6 Further comments on the dual field theory 

As we have remarked above, the string background we have constructed is dual to a flavored 
Klebanov-Witten theory which has gauge group SU(N C ) x SU(N C ), a Landau pole in the 
UV and is conformal in the IR. This is reflected by the behavior of the dilaton e^, which, 
provided the standard orbifold dictionary can be applied to the present setup, is mapped to 
the field theory gauge couplings by 



Air An 

2 ' 2 

9x 92 



2+^2=9^. (5.17) 



Without loss of generality here we consider the case g\ — g\ = g FT = 8ng s e^. The dilaton is 
running and going to infinity for p — > 0, while it is constant in the region p < p q . The beta 
function for gp? can be inferred, holographically, from the dilaton equation 

±e-+ = - 3 -°4M , (5.18) 

dp 47r y 

once the precise radius/energy relation is given. In the massless case, the relation p = 
log (p:/Auv) was used to get the field theory beta function [9]. In our case this relation 
should still be valid with a good approximation when the flavors are very light. Notice that 
this relation suggests that 

p q = \og-^, (5.19) 

A UV 

which, looking at eq. (j5.14p . is just the leading term in the small mass limit where p q —* — oo. 
Within these approximations, we can write the field theory running coupling as 



8vr N 



I 



9ft 47r 



31og- 



A 



uv 




(5.20) 



In this expression, the first term reproduces the one-loop result, while the rest is a string 
predictionjj The power-like factors in the above expression could be interpreted in field 
theory as fractional instantons, or renormalon corrections; the interesting feature is the 
presence of just two such terms. 



6 One has to keep in mind that the string theory calculation of Nf(p) is reflecting a particular choice of 
renormalization scheme in the field theory side. 



6 The Heaviside approximation 

A simple way to model the integrating in/out of the massive flavor degrees of freedom in 
field theory is to forget about the details of the embedding and model the string dual with 
a Heaviside step function Nf(p) = NfQ(p — p q ). Of course this approximation is a source of 
systematic errors, due to the differences between the step function and the true value of Nf(p) 
extracted from the smearing procedure. In the present setup we can measure these errors, 
comparing (see the following sections) physical observables deduced using both the true and 
the approximate supergravity solutions. Estimating these errors is of notable importance to 
have an idea of how well the step function approximation can work in general models like 
[TJ, where the correct Nf(p) may be hard to compute. 

We now present the solution using the Heaviside form for Nj(p), using capital symbols 
for integration constants within this framework. At p < p g , we still have (15.61) . whereas for 
p > p q the solution reads: 

4vr 



e 9 = C 3 ((1 - 6p)e 6 " + C 



. 6 



e f = C 3 ^-6pe 3p [(l-6p)e Gp + C 1/ 
h = 27irg s N c a' 2 (C 2 + f° e'^^dp*) . (6.1) 

As in section 15.31 continuity at p q gives some relations among the different constants: 

br C x = -GpA - e 6 ^ , C, = C 3 (-6 Pq ) l e , (6.2) 



3g s N f p q 

while if we follow again the prescriptions of section IB~4l we get: 

C 1 = C 2 = 0, C 3 = (1 - e 6 P«)- l e . (6.3) 

Under this prescription, the constituent quark mass (15. 5B reads: 



1 8vr e^ 

m„ = J i r ■ 6.4 

q 2Tra'\jg s N f (_ 6 p,)^(l - e 6 ^)i 

Notice that for p q — >• — oo the solutions obtained with the actual Nf(p) approach those 
written above: this suggests that the Heaviside approximation is quite good in the small m q 
limit. 
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7 Wilson loops in backreacted D3D7 models 

In the previous sections we have constructed a string dual of a flavored version of the 
Klebanov-Witten theory with a large number Nf of dynamical flavors of mass m q . Now 
we want to study how these sea flavors affect the non perturbative dynamics of the gauge 
theory. We are going to probe the latter with an external quark- antiquark Q, Q pair with 
mass Mq 3> m q . Due to the presence of a Landau pole in the gauge theory, Mq cannot 
be taken to be infinite and it must lie below the pole. The idea is to study how the static 
quark- antiquark potential depends on the sea quark parameters m q , Nf. The effects of Nf 
flavors on the energy/spin relations for mesonic-like bound states in a similar but localized 
D3D7 system have been studied in [27]. 

The QQ bound state is dual to an open string with the extrema lying on a probe D7-brane 
embedded in such a way that it reaches a minimal distance pq ^> p q from the bottom of the 
space (Mq = 7^-7 !-%, e^ + ^ dp). The string, in turn, bends in the bulk and reaches a minimal 
radial position po- The Minkowski separation L between the test quarks, as well as the total 
energy of the system, depends on p . From this we can deduce the V(L) relation, where V(L) 
is the QQ potential, i.e. the total energy to which the total contribution from the quark mass 
2Mq has been subtracted. The open string embedding is chosen as t = r,y = a, p = p(y) 
where y € [-L/2, L/2] is one of the spatial Minkowski directions. The string worldsheet 
action reads 



s = ~9w / dtd v\l 9tt{p)[g yy {p) + (d y p) 2 g PP (p)] 



2na' 

where g^ v refers here to the string frame metric, which is obtained multiplying (15. ip by ea. 
Defining 

F = y/9tt9yy = e% h' 1/2 , G = y/g t t9pp = e*e f , (7.1) 

we can write [2H EH] the string length and renormalized energy as (the subindex means 
that the quantity is evaluated at p = po) 

fpQ GF 
L{p Q ) = 2 / dp , 

- (PQ Gi 
2vra' ' 

We can use the background solutions found in the previous sections to study how the external 
quark interaction depends on the dynamical massive flavors. In doing so, we vary one of 
the physical parameters Nf,m q ,MQ while keeping the other two fixed. In the numerical 
computations we will present, lengths will be measured in units of y/27Trg s N c a' 2 whereas 
energies in units of (27ro/) -1 (the r that was set to 1 in (13.61) should be reinserted to get the 
correct dimensions). 
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Figure 3: From left to right, plots of V(L) at different values of: g s Nf = 1, 0.6, 0.2, top to 
bottom (at m q = 0.1, Mq = 3); m q = 0.5, 1, 1.5, top to bottom (at g s Nf = 1, Mq = 2.5); 
Mq = 0.6, 1.2, 3, top to bottom (at g s Nf = 1, m q = 0.1). 

From figure [3] we can see that the static quark-antiquark potential is always negative and 
has the same qualitative Coulomb-like behavior as in the unflavored KW conformal case. 
This behavior is not unexpected: in the far IR (i.e. at large L) the quarks are integrated out 
and the theory looks like the KW one. Hence at large L we expect the Coulomb-like behavior 
typical of a conformal theory. Figure [3] also tell us that the absolute value of the potential 
|V(L)| at fixed L is a decreasing function of Nf and an increasing function of m q . This is 
an effect of the screening of the color charges due to the dynamical flavors: the more the 
theory is "unquenched" (large Nf, small m q ) the more the modulus of the quark-antiquark 
force is reduced. Finally, the large L behavior of V(L) is not strongly affected by the choice 
of the cutoff Mq; in the small L region, instead, we see that | V(L)| is an increasing function 
of Mq. We will see in section [7^1 how these results nicely fit with the expected qualitative 
behavior of V(L). 

For all the plots of this section, we have adopted the prescription for integration constants 
given in section \5M It can also be instructive to study how the behavior of the QQ interaction 
is affected by the choice of these integration constants. A brief analysis of the dependence 
on c 2 can be found in Appendix [B] 



7.1 The screening length 

The connected QQ configuration we have studied above is expected to be unstable due to the 
presence of the dynamical flavors. If the QQ string has enough energy, in fact, a dynamical 
qq pair can be popped out from the vacuum with a consequent decay of the heavy QQ state 
into a pair of heavy- light mesonic-like bound states Qq + qQ. The length at which the decay 
can happen is the so called "screening length" L s . In the model at hand, due to the smearing 
procedure, the lighter heavy-light "mesons" are nearly massless. This is due to the fact that 
there exist dynamical flavor branes which intersect the probe brane corresponding to the 
test quarks. The open string stretching between the probe and those dynamical branes has 
nearly massless modes at the intersection. In fact their mass is expected [30] to scale as 
Mq/\/\ (where A denotes the bare 't Hooft coupling), which is parametrically smaller than 
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Figure 4: From left to right, the screening length L s as a function of: g s Nf for m q = 
0.01,0.05,0.1, bottom to top (M Q = 2); m q for g s N f = 1,1.5,2, top to bottom (M Q = 2); 
M Q for g a N f = l,m g = 0.1 (thick), g s N f = l,m q = 0.101 (dashed), g s N f = 1.01, m, = 0.1 
(solid). 



Mq since A ^> 1 in the supergravity regime. The decay into a pair of mesons of this kind can 
happen when the energy of the static configuration, 2Mq + V(L), equals their total mass, 
2M Q /V\. Hence 

V(L s ) = -2M Q (l--L). (7.2) 

Due to the presence of the Landau pole, Mq cannot be as large as we want and it can be 
interesting to study how L s varies with the flavor parameters. Relevant plots (with \/A = 50) 
can be found in figure HI As expected, the screening length decreases as Nf increases and it is 
an increasing function of m q . From the Coulomb-like behavior of V(L) we expect, moreover, 
that L s is decreasing with Mq. The fact that the potential also depends on the cutoff Mq (at 
fixed L, \V(L)\ is an increasing function of Mq as we have seen) does not alter this behavior, 
as it is evident in the figure. 

A recent proposal for improving the quenched approximation in QCD so to take into 
account the screening effects of dynamical flavors can be found in 



7.2 The string breaking length 

It is important to notice that due to the smearing, the decay of the static configuration 
into a given pair of heavy-light bound states is suppressed by l/N c and not by Nf/N c as 
it would happen in a setup with localized parallel flavor branes. A decay rate which is 
only Nf/N c suppressed can be obtained by considering the possibility of decaying into an 
arbitrary sizable fraction of the Nf types of heavy-light bound states. We call the minimal 
QQ separation at which this kind of decay can happen, the "string-breaking length" L sb . We 
choose to define L s & as the length at which a dynamical quark- ant iquark pair with the same 
internal charges as the test quarks can be popped out from the vacuum. The corresponding 
produced heavy-light bound states have a mass given by Mq — m q and thus a large binding 
energy Ej, = —2m q independent on the 't Hooft coupling. The open string describing these 
"mesons" is just stretched along the radial direction of the geometry ( 15. 1ft between the probe 
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Figure 5: From left to right, the string breaking length L sb as a function of: g s Nf, for 
0.03,0.05,0.1, top to bottom (Mq = 2); m q , for g s N f = 0.5, 1,2, top to bottom (Mq = 2); 
M Q for g s N f = l,m q = 0.1 (thick), g s N f = 1.01, m, = 0.1 (solid), g s N } = l,m q = 0.101 
(dashed). 



brane at pq and a parallel dynamical flavor brane at p q [32] . Since the interactions between 
the two produced "mesons" is negligible, it follows from our definition that L sb is just the 
solution of 

V(L sb ) + 2M Q = 2(M Q - m q ) => V(L sb ) = -2m q . (7.3) 

Figure [5] shows that L sb is a decreasing function of both Nf and m q . The latter behavior is 
due to the shift in the meson mass as we vary m q . L sb is an increasing function of Mq which 
tends to flatten at large Mq. 



7.3 Comparison with the Heaviside approximation 

It is instructive to compare the behavior of the Wilson loop observables as obtained using 
the correct string dual solution with g s Nf(p) given in (14.181) . with those obtained using the 
Heaviside step function approximation. In figure [6] we compare the plots for V(L) at different 
values of the parameters. The figure confirms that the step approximation, in the flavored 

V 




Figure 6: The potential V(L) in the step approximation (dashed lines) and with the correct 
Nf(p) function (solid lines), for Mq = 2.5, g s Nf = 1 and m q = 0.3, 1.5 (top to bottom). 

KW model, works quite well when m q is small, or, better said, when Mq — m q is large. 
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Figure 7: Fits of the numerical data (dots) with an analytic formula like (\7.7\i ; on the left, 
m g = 0.01, M Q = 2; on the right, g s N f = 1, M Q = 2. 



Moreover, it is possible to show that the approximation is even better if the step is placed 
where the function Nf(p) has reached half of its value Nf(p step ) = Nf/2. 



7.4 Comparison with field theory expectations 

The behavior of the potential V(L) as well as that of the critical lengths L s and L s b, as a 
function of the flavor parameters, fits very well with a simple analytical result suggested by 
field theory and AdS/CFT arguments. 

We know that in the low energy regime, when the massive flavors are integrated out, our 
model has vanishing beta function. In the string dual description this is reflected into the 
constancy of the dilaton = (fi IR in the p < p q region. The addition of flavors gives a negative 
beta function, whose one loop coefficient behaves as buv ~ ~~ Aff. In the perturbative regime 
m q <C A uv (where A uv indicates the position of the Landau pole), we can match the UV 
and IR couplings at the scale m q so to get the IR value of the coupling 



8tt 2 A uv 



9ir 



m„ 



(7.4) 



This expression is nothing but the field theory rewriting of the stringy result we have found 



for the dilaton in the step function approximation, e^ 



-4ir/(3g s NfPq). Indeed, when 



p q — ► — oo we can safely identify p q ~ \og{m q / A uv ). Now, remembering the expression for 
the QQ potential of a conformal theory with an AdS dual [28J, we expect that for our model, 
in the strong coupling regime 



Vqo(L) 



9 2 irN c 



tNl 



1 



(7.5) 



in the large L limit. From this expression we can deduce, in the same limits, the screening 
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length L s 



L s « l ^ Wc = , (7.6) 



and the string breaking length L^ 



'sb 



m q JN f log 



(T.7) 



A 



j_i_ 



The above analytical functions share with our numerical results the same qualitative behavior 
with the flavor parameters. The expression for the screening length is not expected to fit 
well with our data, since the approximation we used to get eq. ( 17.61) requires L to be large. 
Instead, in figure [7] we have fitted our numerical data with a formula for L sb as in eq. (17.71) 
- treating Aw and an overall constant as parameters - finding a very good agreement. 
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A D7 embed dings from the orbifold 

There are two main classes of holomorphic embeddings for D7 branes on the conifold: class 
I has Z\ — ji [22] as representative embedding equation, while class II is represented by 
z\ — z 2 = fi [23]. The two classes are argued to correspond to two different ways of adding 
flavors to the KW model. Let us focus here on the massless \i = case. Class I embeddings 
have two branches of D7 branes, each branch adding fundamental matter to the first node of 
the KW quiver and antifundamental matter to the second. In the class II case a D7-brane 
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provides fundamental and antifundamental matter to one node. Since the KW theory can be 
obtained starting from the J\f = 2 conformal theory describing N c D3-branes on the C /Z2 
orbifold singularity [10] it could be useful to classify the possible kinds of D7 embeddings 
starting from the orbifold picture. This could also help in finding the correct superpotential 
terms corresponding to the chosen embeddings. A discussion on these issues can be found 
in [22]. 

The orbifold theory is a quiver with gauge group SU(N C ) x SU(N C ) and, in terms of N = 1 
components, the same bifundamental multiplets A^ B t as the K W theory, plus two adjoint 
supermultiplets $1 and $2, one for each of the gauge groups. The superpotential of the 
theory schematically reads 

W orb = $i(AiBi - A 2 B 2 ) - $ 2 (BiAi - B 2 A 2 ) . (A.l) 

There are mainly two ways of adding massless flavors to both nodes of the quiver: either we 
add to the above superpotential cubic terms coupling the new flavors to the bifundamentals, 
i.e. (modulo SU{2) x SU{2) rotations) 

W orbI = W orb + q 1 (A 1 - A 2 )q 2 + q 2 (B 2 - B x )q x , (A.2) 

or we couple the new flavors with the adjoint supermultiplets 

Worbii = W orb + qi$iqi - q 2 $ 2 q 2 . (A.3) 



These two possibilities have been considered in the T-dual Type IIA picture [331 El] to which 
we refer for the sign conventions. The superpotential in (IA.2I) reduces supersymmetry by one 
half, while that in ( 1A.3I) does not break the M = 2 supersymmetry (the corresponding string 



model is just an orbifold of the D3D7 system in flat space). In case I each branch of D7 branes 
extends only along two directions inside the orbifold. In case II, the D7-branes are taken to 
be extended along the 4 orbifolded directions. It is easy, in this case, to see that the flavored 
theory has a Landau pole: the one loop coefficient of the beta function (which does not get 
perturbative corrections beyond one-loop due to the M = 2 supersymmetry) for each group 
is equal to b = 3N C — N c — 2N C — Nf = —Nf (where the negative terms — N c , —2N C , —Nf 
are the contributions of the adjoint, bifundamental and fundamental matter superfields, 
respectively) and on each node the theory is thus IR free. 

Let us now go to the J\f — 1 conifold theory by adding to the above superpotentials the 
appropriate mass term for the adjoints 

W m = ^(*?-*2). (A.4) 

Integrating out the adjoint fields in case I, will produce a superpotential of the form 

Wi = W KW + h qi(A t - A 2 )q 2 + h 2 q 2 (B 2 - B 1 )q 1 , (A.5) 
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while in case II one will get 

Wn = W K W + k q\ [AiB l - A 2 B 2 ] qi + h 2 q 2 [ B i A i ~ B 2 A 2 ]q 2 + k x (q iqi ) 2 + k 2 (q 2 q 2 ) 2 . (A.6) 

Writing the effective mass terms for the flavors in a compact matrix form qMq, the equation 
det M = gives, for the case I, A i B 1 + A 2 B 2 — A X B 2 — A 2 B X = 0, i.e. z\ + z 2 — z 3 — z& — 0. 
This is nothing but one of the possible rotations of the z\ = embedding: in fact it is our 
generalized embedding (13.111) with a 3 = a 4 = — 1. In case II, det M = gives A X B X — A 2 B 2 = 

0, i.e. the z\ — z 2 = embedding equation. 

The massive generalizations of the embeddings above correspond to the addition of stan- 
dard mass terms miqiqi + rn 2 q 2 q 2 to the massless superpotential. In case I the mass terms 
break the classical flavor symmetry group to U{Nf), while in case II, the U(Nf) x U(Nf) 
symmetry is preserved. 

It is also tempting to write the class I D7 embeddings using the orbifold coordinates. The 
C 2 /Z 2 orbifold is described by the xy = z 2 equation in C 3 . The massless embeddings of class 

1, should thus be described by an equation like z — 0, giving rise to the two branches x = 
and y = corresponding to D7 branes extended only over one half of the orbifold directions. 

B Dependence on the integration constants 

In Section[5]we have constructed the background solution obtained with a precise prescription 
for the integration constants Ci, 02,03. We would like to briefly discuss here what would 
happen with a different choice of constants. As explained in Section 15.41 c\ 7^ would 
produce a singularity in the IR which we want to avoid. On the other hand, varying C3 just 
produces a rescaling of e* , e 9 , h and, for the static quark-antiquark system, a rescaling of L 
and V which does not change the qualitative form of the V(L) curves. Finally, c 2 < would 
make the geometry singular at some point below the Landau pole, a behavior that we also 
want to avoid. Let us thus focus on modifying our original prescription c 2 = to c 2 > 0. 
The behavior of the V(L) static potential is crucially affected by the value of c 2 , as shown 
in figure For every c 2 > the large L (p — > — 00) regime of V(L) is unchanged: the 
potential approaches zero from below since the theory is IR conformal. However, above a 
certain critical value of c 2 (which depends on the other physical parameters), one observes a 
double turnaround of the V(L) graph at shorter distances: the potential has a discontinuity 
in the first derivative, and there is a first order phase transition. 

This kind of behavior was first found in [35] and afterwards discussed in several different 
contexts [361 Ej- I n ^ ne region of L where V(L) is triple- valued, the point with lowest V is 
stable, the intermediate one is metastable and the one with larger V(L) is perturbatively 
unstable (always referring to longitudinal perturbations of the string). Stability issues of 
this kind of strings were discussed in detail in [37]. 
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Figure 8: Plots of V(L) with c 2 = 10, m q = 0.2, g s N f = 1, Mq = 2.5. The dashed (resp. 
solid) line refers to the result obtained using the Heaviside approximation (resp. the correct 
solution with Nf(p)). 

The double turnaround region is visible also in the massless m q = case. As it is evi- 
dent from figure [HI the Heaviside approximation does not modify qualitatively the double 
turnaround behavior of the potential. 
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